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Abstract. For the spherical symmetric case, all quantities describing the relaxed dark matter halo can be expressed 
as functions of the gravitational potential "1>. Decomposing the radial velocity dispersion ar with respect to $ at 
very large and very small radial distances the possible asymptotic behavior for the density and velocity profiles 
can be obtained. If reasonable boundary conditions are posed such as a finite halo mass and force-free halo 
centre the asymptotic density profiles at the centre should be much less steep than the profiles obtained within 
numerical simulations. In particular cases profiles like Plummer's model are obtained. The reasons of that seeming 
discrepancy with respect to the results of N-body simulations are discussed. 
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1. Introduction 

The formation of structures by gravitational interact- 
ing cold dark matter (CDM) is one of the outstanding 
paradigms in cosmology. The luminous baryonic matter is 
embedded in dark matter halos. Numerical studies of the 
structure formation, which allows only for gravitational 
interactions, predict the distribution of galaxies and clus- 
ters of galaxies in excellent agreement with observations. 
Special attention has been devoted to the investigation 
of halos as typical aggregates forming after collapse and 
violent relaxation. 

Early studies seemed to indicate that not only the 
distribution of halos but also their densit y profiles de- 

uinn et al.1. 
199lt 

However, numerous numerical inves- 
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tigations have shown that halos have an almost univer- 
sal profile that is independent of the mass of the halo 
ijNavarro et all ri996. 1997), the initial fluctuation spec- 
trum and the underlying cosmological model. Although 
the resolution of the simulations has increased dramat- 
ically since the first studies, high resolution simulations 
are still not able to provide a reliable de termination of 
the s hape of the innermost density profiles. iNavarro et alJ 
l)l99& |1997 ) first proposed an analytic approximation 
for the density profiles obtained in simulations. They 
obtained p oc (r/rs)~^(l -I- r/rs)~'^, where denotes 
the scale radius of a given halo. This NFW-profile im- 
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plies an asymptotic behavior p oc r^" with a power 
index of n = 1 for the density profile throughout the 
central halo region. Other numerical studies produced a 
signi ficantly steeper inner slope corresponding; to n ~ 
1.5, ('Moore et al.. 'l998'. '1999*: Ghigna et all 'l998t Eoocj 
j:ukushige & Makino, 2001 : Klvpin e t al.. ,2001') . In con- 
trast to these results IPower et al.l l)20o'^ did not find 
an asymptotic power behavior with fixed index at all 
but rather a continuously decreasing index with values 
down to 1.2 at the innerm ost radius they could resolve. 
Finally, 'jing fc Sutol l)2000() claimed that the asymptotic 
behavior of the innermost region depends on the mass of 
the halo. 

The observations yield a different picture: From the 
rotation curves of stars the density profile of the un- 
derlying dark matter halo of a galaxy can be inferred. 
The observation that Low Surface Brightness (LSB) and 
High Surface Brightness (HSB) galaxies follow the same 
TuUy-Fisher relation requires (in the conventional pic- 
ture) that LSB galaxies arc domin ated by dark r natter 
(DM), e.g. .Verheiien ( 19971. de Blok fc McCTaugh I lll99fJ) 
. They are promising candidates for determining the in- 
nermost density profile. A 'cuspy', n > 1, dark matter 
core is ruled out by many studies ijMcGaugh fc de Bloki 
ll998Hde Blok et aLlEool . It is also rule d out by modeling 
the in ner rotation curve of the Galaxy ijBinnev fc Evana 
l200lh . In contrast to this results some recent inves- 
tigations of dwarf galaxies suggest that an innermost 
slope of n » 1 is con sistent with rotation curve data 
l)van den Bosch fc Swaters, .200^ . Thus, present observa- 



2 



J. P. Miicket and M. Hoeft: Density Profile Asymptotes at the Centre of Dark Matter Irtalos 



tions p ermit an inner slope of n < 1 at most. lEl-Zant et al.l 
l)200l|) suggest that a 'cuspy' halo core may be eliminated 
by dynamical friction of an initially clumpy gas distribu- 
tion. This indicates that the gas content of a halo may 
aff ect the dark matter profil e. 

iTavlor fc Navarrol l)200l[l argued that the recurrent 
merging process results in a phase-space density pro- 
file which decreases according to a power-law. With 
pjci^ oc y-i- STS ^];^gy reprod uced well the numerical den- 
sity profiles. iNavarrol l|200lf) argued that the entropy in 
a halo is most uniform distributed under the condition 
of a power-law phase-space density. This results in ob- 
served density profiles. In contrast to an universal den- 
sity profile some authors predict again a dependency 
on the for mation history a nd the underlying cosmolog- 
ical model llSver fc Whitl Il998t iNiisser fc Shethl. Il999t 
iLokas fc HoffmanL 1200(11)7 Recen t p apers bv iDekel et al.l 
1120021) . IWeinberg fc Katzl ll2002t) and lSellwoodI ll2002l) dis- 
cuss the influence of subsequent mergers or central bars 
onto the formation of cusps at the dark matter halo cen- 
tre. Though the discussion about the bar influence is still 
controversial the common outcome of those investigations 
is that a cusp seems not to be a generic feature of a re- 
laxed dark matter halo but is rather the result of further 
interaction with the halo surroundings. The effect on the 
halo profiles due to the interaction between the gas and 
the dark matter on the profiles is still an open question. 
The possible structure of pure dark matter halos is still 
a basic problem and needs to be considered separately to 
get insight into the more general behavior of the combined 
system of dark matter and gas. 

In this paper we suppose that the halos have spherical 
symmetry and are sufficiently well relaxed, i.e., the halos 
can be described by the steady state Jeans equation. We 
systematically investigate the asymptotic behavior of the 
quantities describing the halo at radii comparable to and 
larger than the virial radius (outer region) and near the 
central region at small radii. The assumptions to be made 
will be kept very general, e.g., demanding the finiteness of 
the halo mass. 

The most essential and probably the most restrictive 
boundary condition we will adopt is the demand for a 
nonsingular potential throughout. This appears as an ap- 
propriate assumption for real dark matter halos if the ex- 
istence of a central point mass is forbidden. 

A particular aim of the present paper is to investigate 
under which conditions a power asymptote with respect 
to the radius exists for the density profiles of dark matter 
halos if approaching the very centre. 

2. The asymptotes at large radii 

The system of equations describing spherically symmetric 
halo configurations consisting of coUisionless dark matter 
particles is formed by the Jeans equation and the Poisson 
equation. The halo structure is completely described by 
the density p(r) , the components of the velocity dispersion 
(7^ = cr^ -I- CTg + a\ and the gravitational potential $(r). 



Due to spherical symmetry it follows = cr^. In order to 
close up the above system of equations one needs two more 
relations. Even if supposing a nearly constant anisotropy 
parameter ^ —\ — CTg/a^ still a further relation is needed. 
Asymptotic considerations may provide such a relation 
from some basic assumptions. 

For the spherical-symmetric case, the potential $(r) 
is an increasing monotonic function with respect to the 
radius r; d$(r)/(ir = GM{r)/r'^ > at r > , where 
M(r) is the cumulative mass at radius r. Therefore, all 
quantities can be expressed as functions of <I> instead of r. 
The potential $ is determined up to an arbitrary constant 
which has to be fixed by the chosen boundary conditions 
usually given either at r — > oo or at r = or at a finite 
boundary r = i? if it exists for the considered object. We 
suppose a non-singular potential throughout. Then it can 
be recalibrated to fulfill <I'(0) = at the center of the halos. 
Further, we assume the potential to be sufficiently smooth 
to fulfill the boundary condition for d^/dr at r = 

d^ 

(1) 

which we adopt throughout. 

We suppose that the entire mass of the halo is finite. 

The latter condition means that at r — > the density 
p must not increase faster than p oc r'^~^ where e > can 
take an arbitrary small value. The condition requires 
that the density p at r — » must not increase faster than 
p{r) oc r^^+* with e > 0. Otherwise d^/dr = M{r)/r'^ 
becomes singular at r — > 0. Thus the condition Q leads 
to an even stronger restriction for the density behavior 
fulfilling the supposition of finite halo mass at r au- 
tomatically. 

At r ^ oo the potential $(r) approaches ^oo- Since 
at r ^ oo the probability to find particles with nonzero 
velocities should vanish, it follows also cr^ ^ at r — > cx). 
A decomposition of tr^ into a power series in the vicinity 
of $oo can be achieved, where the leading term is oc ($oo — 
with p > 0. 

Therefore, af. can be represented in terms 
al = a($oo - W = a-^^ir). (2) 

where a is a positive constant which has to be determined. 
Here we adopted th e so-ca lled relative potential (see e.g. 
iBinnev fc Tremaind (|l983)) * = «'oo - $(r). 
In this case the Jeans equation 

— h 2/^-0-^ -p—. (3j 

dr r dr 

can be integrated and we obtain 

p oc ^If-Pr-^f^ exp( J (4) 

a(l -p) 

In addition to a finite halo mass the density must van- 
ish p — > at r — > oo. Although the asymptotic behavior 
at '5 is different for p < 1 from the behavior for p > 1 
this is satisfied by all possible p. The singular case p = 1 
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must be considered separately. If p < 1 the asymptotic 
behavior of p{r) at ^ is determined by the terms in 
front of the exponential term 



Inserting © into the Poisson equation 

- — (r2 — ) = 47rGp, 
ar ar 



we get 



1 d 



dr dr 



) + AttGc-^' 



-2/3 _ 



= 0, 



(5) 



(6) 



(7) 



One solution of Eq. Q is a power law with respect to 
r which gives at the same time the asymptote for large r 
where the density p must be a decreasing function. 

Assuming ^ oc r' a relation between the exponent 
I and the parameters /3 and p introduced above can be 
given: 

2(1-/?) 



I 



(8) 



Since always /? < 1 it follows from that Eq. {Tj) does 
not permit a power asymptote which fulfills the condition 
^'(r) ^ at r ^ oo. 

For the cases p > 1 according to Q the density 
goes exponentially fast to zero for 5" — > Q at r oo . 
Those models were obtained, e.g., by iMilgromI ((l996![). 
This means that the density term in the Poisson equation 
vanishes much faster than the first two terms d^'^/dr^ + 
(2/r)d'^ /dr. Hence the asymptotic solution for 5* is de- 
termined just by these terms and we obtain ^'(r) oc l/r. 
The solution with exponential decrease of the density 
at r ^ oo describes a halo the mass of which is dis- 
tributed within an finite radius i?eff and the asymptotic 
solution for is equivalent to the potential describing 
the outer region at r > i?eff approaching the outer so- 
lution for a spherical symmetric finite body. The radius 
i?eff is determined by the potential value when the ex- 
ponent in becomes effective. Provided the halo mass 
Mq is known then the halo radius can be estimated to be 
i?eff «GMo(a(p-l))i/(f-i). 

Now, let us consider the particular case p = 1 . In this 
case, the Jeans equation can be integrated neglecting ini- 
tially a possible anisotropy of the velocity dispersion. We 
find 

P = c^\ (9) 

where 7 (1 — a) /a. We use Eq. to eliminate p 
from Poisson's equation and find 



dr dr 



(10) 



which is known as Lane-Emden equation. One solution 
of Eq. IjlOf) is a power law with respect to r 



Assuming p oc r", a relation between the exponent n 
and the coefficient a can be obtained: 



2(1-^ 
2a - 1 

Since a > and n, < we get either 
< a < ^ 



or 



a > 1 



(12) 



(13) 



(14) 



The latter condition can be excluded since it leads to in- 

< i satisfies the 
00, at the same 



finite halo masses. The condition < a < ^ satisfies the 



boundary condition ^'(r) — > at r 
time. 

The demand of a finite halo mass will be satisfied if 
71 < —3. This restricts a further via 



1 1 

T < " < 7: 
4 2 



(15) 



The special case of a = 1/4 results in the asymptotic 
behavior of the NEW profile. 

If we admit anisotropy for the velocity dispersion 
which is indicated by the results of the halo analysis from 
N-body simulations for the outer region then we have to 
consider the Jeans equation with /3 ^ 0. 

Assuming (3 ~ const, and using the relation (j^J we are 
again able to integrate ((SJ obtaining 



\I/7 



(16) 



The asymptotic power laws for 5* and p are then 5* oc 
r" with n — 2(1 — /3)/(l — 7) and p oc r™ with m = 
nj — 2(3 = 2(7 — /3)/(l — 7). Arguments as given above 
lead to the condition for a 



1 1 
< a < - 

-2(3 2 



(17) 



(11) 



For /3 = 0, the above obtained results are recovered. 
The case (3=1 must be considered separately. We solve 
the Jeans equation for this case and obtain that a = 1/2 
is necessary for the existence of a power asymptote. The 
finite-mass condition rules out this asymptote, however. 

The case a = 1/2 is also not included within the above 
considerations. For an arbitrary /3 7^ 1 this leads to the 
asymptotes of Bessel functions at r ^ 00 also not fulfilling 
the finite-mass-condition. 

Thus, an anisotropy with respect to the velocity dis- 
persion (3^0 confines the possible range for a, i.e., this 
leads to the more restrictive condition H17() . 

3. The asymptotes at the innermost halo region 

We shall now consider the central part of the halos in 
more detail. The knowledge about this region is uncer- 
tain partly because of the lack of resolution and also due 
to the poor particle statistics when approaching r 0. 
So the question whether the solution approaches a power 
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1 10 100 1000 



Radius R (in arbitrary units) 

Fig. 1. The figure shows the shapes of the obtained 
asymptotes at r ^ oo covering the area between the 
asymptote with {p = l,a — > 1/(4 — 2(3)) (thick solid 
line) and the dashed vertical line (p = l,a 1/2). For 
comparison two asymptotes for p > 1 are shown: The 
asymptote with p = 1.3 approaches for decreasing p the 
[p — l,a — » 1/(4 — 2/3) asymptote (thick dashed hne), 
the p = 3 asymptote exhibits an extremely steep decrease 
(dotted line) . Asymptotes with increasing p approach the 
vertical dashed line. For p > 1 the density falls off very 
rapidly after some characteristic radius and describes a 
finite halo. All asymptotes are given for (3 — 0.8. 

asymptote within the central region or not will be consid- 
ered. Provided a power asymptote exists then its possible 
properties will be analyzed. 

As was argued above for the spherical-symmetric case 
all quantities could be represented as functions of the po- 
tential <i>(r). For the above adopted boundary conditions 
($(0) = 0) we consider the decomposition of the velocity 
dispersion into a power series with respect to <i>. In order 
to allow for the case that (7^ exhibits a power-like behavior 
with an arbitrary exponent at <1> = we consider the more 
general decomposition 

^2 ^ Co$''(l + ai$ + a23>^...) (18) 

We first consider the case k 7^ 0. Then the leading term 
at r is 

al = Co$'' (19) 

where k > if supposing regularity of a'^(r — > 0). 
Neglecting anisotropy, we get from the Jeans equation: 

poc$-'-exp(- " ) (20) 
Go(l - k) 

In order to avoid an exponential-like singularity of the 
density at r ^ for k must be valid k < 1. The case k > 1 
would also lead to an infinite halo mass as long $ can be 
represented as some power of r at r — > 0. Although we will 



restrict ourselves on just this case in some considerations 
below, we should be are aware of this is an additional as- 
sumption. Therefore, the density behavior near the center 
is given by 

(21) 

where $ — > 0. 

This is true as long as the exponent in 120|) is small 
with respect to unity, i.e., $ < (Co(l - k))^/^^"''). The 
potential $ can be scaled with respect to its value $00 at 
r ^ 00. Then Cq = c^"^" where the comparison with 
the result s of high-re s olutio n numerical N-body simula- 
tions fsee lHoeft et al.l l)2003l) l provides a numerical value 
of about c ~ 0.3. 

Inserting Eq. H21(l into the Poisson equation we can 
find the asymptotic solutions for $ at r — > 0. Due to non- 
linearity several kinds of solutions are possible, in prin- 
ciple. If we restrict ourselves to solutions for $ having a 
power-law behavior with respect to r near r = then we 
get for the asymptotic behavior 

$ oc r2/(i+«) (22) 

cr^ oc cx r2«/(i+«) (23) 

Though p cx <^^'^ cx 7.-2k/(i+k) jg singular at r ^ 
the halo mass remains always finite if < k < 1. 

The case k = 1 must be considered as a special case 
solving eq. fH^ and leads to a separate solution not covered 
by Ian 

p oc <i>-(^) oc r-2(i+Co)/(2Co+i), (24) 

From ^ follows the condition 

2(l + Co)/(2Co + l) < 1 (25) 

which cannot be satisfied for any possible Cq. Therefore 
the case k = 1 is ruled out. 

An almost isotropic velocity dispersion is justified at 
r — > by the results of numerical simulations. In general 
/9 is a function of r, which is small or vanishes if r 0. 
From the simulation results, it is known that (3{r) is an 
increasing function of the radial coordinate. If at r — + /? 
is non-zero, i.e., /3 ^ /3o at r — > 0, then 

p(xr-2^"<i>-'-exp(- f ' " ) (26) 
Co(l - k) 

$ oc r2(i-/5o)/(i+«) (27) 

<jI oc oc (28) 

For /3 — > at r — > we can make a power law ansatz 
P oc with p > 0. This leads to a solution for p similar to 
the expression (|20|l containing in the exponent an additive 
term proportional to — r'^ which leads at r — > again to 
the case (|21ll . If we pose the additional condition that the 
term d{pa^)/dr being the analog for the pressure gradient 
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0.010 0.100 1.000 10.000 

Radius R (arbitrary units) 

Fig. 2. The figure shows the range of possible density 
asymptotes at r — > (upper left shaded area) . This area is 
bordered by the limit case asymptote with (k = 1/3, m = 
1/2) and the asymptote {k — 0,m = 0). The dashed dot- 
ted vertical line denotes the radii up to which the approx- 
imation p oc reliably works. For larger radii up to the 
right shaded area the expression (|20|l must be taken into 
account. The thick solid and the dotted lines give the cor- 
responding behavior. The thick dot-dashed line gives the 
behavior of the effective power index (the absolute value 
is plotted). Throughout the range of radii where $ can 
still be considered as small the index |m| is increasing and 
of order unity. 



considerations lead to a restriction for the introduced pa- 
rameter e, namely it must be e > 1/2, which exactly cor- 
responds to the above obtained relation with respect to 
m. 

The expression H2()|l for p indicates that at radii signif- 
icantly apart from zero (^/Cq^*-^ is not longer small) 
the density falls off more steeply than the obtained power 
asymptotes (see Fig|2l . This may explain that the asymp- 
totic exponent m = 1/2 is not as steep as the exponents 
found in TV-body simulations. The latter ones can be ob- 
tained only at radii still sufficiently far away from ?' = 0. 

We have to stress, that the obtained asymptote H22|l 
at r — > is the solution if supposing a power-law behavior 
for $ with respect to r. If replacing the density p in the 
Poisson equation by its expression (|20|l we obtain a highly 
nonlinear differential equation which posses also solutions 
leading to a non-singular behavior of the density at r — > 
as given below. 

The result of the above considerations can be expressed 
as follows; If one supposes a behavior of the quantity p 
obeying a power law with respect to the radial coordinate 
p oc r~™ then the power index is restricted hy < m < 
1/2 at r « 0. 

Note however, that even if we are not demanding that 
<i> is a power law with respect to r the more general relation 
(O stiU holds. 

We now consider the case k = which represents the 
possibility of a non- vanishing velocity dispersion at r = 0. 
Taking into account in the Jeans equation for only the 
linear contributions of $ and d^/dr we obtain for the 
density 



should vanish or should be regular, at least, at r — > then 
it follows immediately (3{r ^ 0) — » 0, i.e., also /3o = 0, 
which is consistent with other considerations. 

We use the relation d{pa'^)/dr = —pd^/dr and sup- 
pose that p is a power law with respect to r and obtain 



p oc (1 -I- ai<I> 



(30) 



d{pa^ 
dr 



„(-3K+l)/(l + fe) 



(29) 



Note, in order to obtain the correct expression being 
consistent with the right-hand side of the Jeans equation 
one has to perform first the derivative d{pa^)/dr using 
the expression (|2()|l and considering the limit $ ^ only 
afterwards. 

Regularity of d{pa'^)/dr at r ^ leads immediately 
to the condition < k < 1/3. Thus the behavior of p 
at r « is given by p oc r~™ with a possible range for 
< m < 1/2. 

This result is a consequence of demanding the right- 
hand side of the Jeans equation to be non-singular at 
r — > 0, i.e., requiring not only d^/dr — but also the 
product of pd^/dr to be finite, at least. This leads in- 
evitably to a further restriction of the density behavior. 
One can obtain the above relation straightforward by do- 
ing analogous considerations as with respect to the finitc- 
ness of halo mass resulting from the condition Q). These 



From the Poisson equation we get then the asymptotes 
for $ oc . According to eq. H3U|) the density remains finite 
at r ^ 0, in this case. The shape of the density asymptote 
in the vicinity of r — > is determined by ai and Cq. E.g., 
for ai = 1/3 and Co = 2 the asymptotic density profile is 
the Plummer's law. 

Note, for all cases where > 1 coefficients in (|18|l 
vanish, i.e., ai = 0,i — I, ...N we are led to the standard 
Emden form and the approximate solution describes an 
isothermal sphere. However the condition $(0) — leads 
inevitably to $(r) oc at r 0. All those solutions lead 
to an asymptotic core with constant density. The larger 
the non- vanishing dispersion ctq at the centre is the faster 
the core solution will be attained. 



4. Summary and conclusions 

We have systematically investigated the asymptotic be- 
havior of the quantities describing spherically symmet- 
ric relaxed dark matter halos. We have drawn partic- 
ular attention to the behavior near the centre r — 0. 
The most restrictive boundary condition we have used is 
$(0) = d^/dr{0) = 0. This means that the dark matter is 
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considered to be continuously distributed, i.e., the gravi- 
tational potential at r = remains finite and no singular 
point mass at the centre is allowed. 

This is in agreement with the suppositions for deriving 
the Jeans equation from the coUisionless Boltzmann equa- 
tion. It is demanded that each particle is moving within 
the smooth potential of all other particles and that no 
close two-body encounters take place. The latter is cer- 
tainly true for the real dark matter distribution. Within N- 
body simulations encounters cannot be entirely excluded. 
The probability for an encounter between two particles is 
proportional to the particle density. Thus especially near 
the dense halo centres infalling particles with small im- 
pact parameter may undergo "collisions" just within this 
region. Therefore, one may generally expect some devia- 
tion of the particle distribution in the halos formed during 
N-body simulations from the results given here where we 
assume that the suppositions for the Jeans equation are 
well satisfied. 

The asymptotic behavior of the velocity dispersion as 
function of the potential $ at small and large radii can be 
given by ((TH)l and respectively. Provided the asymp- 
totic density profile at large r is characterized by a power 
asymptote with respect to r then the parameters enter- 
ing the asymptotic relation between the velocity disper- 
sion and the gravitational potential cr^ = Q;(<i>oo ^ 'i')^ can 
be further restricted. For p < 1 no power law asymptote 
does exist which fulfills the boundary condition p — > 
at r — > oo. For p > 1 the density profile exhibits an ex- 
ponentially fast decrease at large distances and the halo 
structure is given by the mass concentration within a finite 
radius. For the case p = \ the parameter a is restricted by 
the relation 1/(4 — 2/3) < a < 1/2. This is in agreement 
with the observed velocity dispersion profiles. 

We get at small radii an asymptotic behavior which 
seems to differ from the behavior derived from current 
results of numerical simulations. Assuming a power law 
behavior of the density profile p(r) oc r^™ at r ^ our 
analysis predicts m < 1 for all cases being in agreement 
with the initial conditions. If we demand that the gradient 
for the pressure analog pc^ should vanish at r = 0, then 
this leads to the even more restrictive relation m < 1/2. 
Thus, if at smaU radii where <i> < (Co(l - k))^/^^^'') is 
valid, a power asymptote p(r) cx r~"^ exists then m must 
obey the above restrictions, i.e., a power asymptote ac- 
cording to the NFW profile cannot be continued to the 
innermost region under the assumption that the halo can 
be described by the Jeans equation. 

While an analytic behavior of with respect to $ 
at $ ^ can be well justified, the existence of a power 
asymptote with respect to r for all the quantities is not 
mandatory. In particular, postulating the existence of a 
power asymptote for the density with respect to r might be 
a too strong requirement. All solutions with nonvanishing 
velocity dispersion at the very centre lead asymptotically 
to a core with p(0) = po < oo. 

We see two reasons that our results seem not to be in 
agreement with the results of the numerical simulations. 



(i) The resolution of the available numerical simula- 
tions is probably not yet high enough to resolve the region 
$« (Co(1-k))1/(i-'=). 

However, as can be seen from the relation H2U|I the ef- 
fective power index becomes steeper with increasing radii. 

This is illustrated by Fig. 12 where the absolute value 
of the effective power index |m| is shown to be increasing 
to values of about and larger than unity although it is still 
^/Cq^*"^ < 1. At $ > 1 the approximation is not longer 
valid. Thus the obtained density asymptote approximates 
the NFW profile at the low-radius-end. This would ex- 
plain why the numerical results obtained for radii where 
<I'^~'^/Co/(l — k) -C 1 is not longer valid lead to a much 
steeper profile (m > 1). 

(ii) An other reason could be that the halos having 
formed within the simulations are not yet sufficiently iso- 
lated from their surroundings. Therefore, a description by 
the steady state Jeans equation may be partly question- 
able. 

Mergers and further matter inflow may alter the phys- 
ical situation essentially and can lead to a steepening of 
the inner density profile IPekel et aTl l)2002() . 

On the other hand a couple of observations seem to 
support strongly the existence of an asymptotic density 
profile at r — > being much less steep than m = 1, e.g., 
dc Blok ct al. ( 2001). 

This may indicate that the inner density profile near 
the centre depends on whether or not the considered halo 
is interacting with the surroundings and on the way this 
happens. 
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